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We show that when a self-avoiding polymer chain is pulled off a sticky surface by force applied 
to the end segment, it undergoes a first-order thermodynamic phase transition albeit without phase 
coexistence. This unusual feature is demonstrated analytically by means of a Grand Canonical 
Ensemble (GCE) description of adsorbed macromolecules as well as by Monte Carlo simulations of 
an off-lattice bead-spring model of a polymer chain. 

Theoretical treatment and computer experiment can be carried out both in the constant-force 
statistical ensembl whereby at fixed pulling force / one measures the mean height (h) of the chain 
end above the adsorbing plane, and in the constant-height ensemble where for a given height h 
one monitors the resulting force (/) applied at the last segment. We find that the force-assisted 
desorption undergoes a first-order dichotomic phase transition whereby phase coexistence between 
adsorbed and desorbed states does not exist. In the /-ensemble the order parameter (the fraction of 
chain contacts with the surface) is characterized by huge fluctuations when the pulling force attains 
a critical value fo- In the ft-ensemble, in contrast, fluctuations are always finite at the critical height 
ho- 

The derived analytical expressions for the probability distributions of the basic structural units of 
an adsorbed polymer, such as loops, trains and tails, in terms of the adhesive potential e and /, or h, 
provide a full description of the polymer structure and behavior upon force-assisted detachment. In 
addition, one finds that the hitherto controversial value of the universal critical adsorption exponent 
(j} depends essentially on the extent of interaction between the loops adsorbed chain so that may 
vary within the limits 0.39 < (f) < 0.59. 

I. INTRODUCTION 

Over the past decade, experimental force spectroscopy techniques such as Atomic Force Microscopy (AFM) and 
optical or magnetic tweezers emerged as novel methods which allow the manipulation of individual polymers with 
spatial resolution in the nm range and force resolution in the pA^ range [llQ- One can thus study the mechanical prop- 
erties and characterize the intermolecular interactions of a single macromolecule which leads to better understanding 
of the material elasticity on a molecular level 3, 4], enables measuring the receptor - ligand binding strength^, or the 
determination of friction- induced energy dissipation during the movement of a macromolecule on a solid surface [gJ. 

The rapid development of experimental techniques has been followed by theoretical considerations, based on the 
mean - field approximation 0], which provide important insight into the mechanism of polymer detachment from 
adhesive surfaces under external pulling force. A comprehensive study by Skvortsov et al. [8!] examines the case of 
a Gaussian polymer chain. One should also note the close analogy between the forced detachment by pulling and 
the unzipping of a double - stranded DNA. Recently, DNA denaturation and unzipping have been treated by Kafri 
et al. |9j] using the Grand Canonical Ensemble (GCE) approach [l^, |Tl| as well as Duplantier's analysis of polymer 
networks of arbitrary topology . An important result concerning the properties of adsorbed macromolecule under 
pulling turns to be the observation that the universal exponents (which govern polymer loops statistics) undergo 
renormalization when excluded volume effects between chain segments are taken into account. In this work we use 
similar methods to describe the structure and detachment of a polymer chain from a sticky substrate under pulling 
and demonstrate the unusual properties of this phase transformation in two conjugated statistical ensembles. 

II. THEORY OF CHAIN DESORPTION 

A. A simplified case of detachment 

In order to illustrate the problem with chain detachment under pulling, we start with a simple example, cf. Fig. [2 
which shows schematically a case when N — m chain monomers are adsorbed on the plane while the remaining m 
monomers form a stretched tail subjected to external force /. Consider for simplicity a phantom chain with no 
excluded volume interactions between the segments. The partition function of such Gaussian chain can be written as 
n{m) — exp[e(A^ — m) — mfa/kBT] where /i2 denotes the so called connective constant in d = 2 dimensions 
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(e.g., = 2.6 on a cubic lattice). The dimensionless adsorption energy e = e/kBT measures the energy gain per 
contact with the surface while the work to detach and move m beads a distance a away from the plane is af /ksT. 

Evidently, the corresponding free energy F/ksT = — hiil{m) w — iV(ln /i2 + e) + m[(ln /i2 + e) — fa/ksT] grows or 
declines with varying to, depending on the sign of the expression in square brackets. Therefore, one can readily define 
a critical detachment force /z?(e) = kBT{\npL2 + e)/a such that for f < fo one finds a minimum of at to = (the 
chain is completely adsorbed) whereas for f > fo the lowest free energy is reached for to = whereby the polymer 
is entirely detached from the surface - Fig. [TJ At the critical value f = fo the free energy becomes independent of to, 
indicating even within this oversimplified consideration (which neglects the presence of loops in the adsorbed state) 
that any number of chain contacts with the adsorbing plane becomes equally probable. Evidently, by just crossing 
the critical line foi^) the polymer chain undergoes an abrupt transition between an adsorbed and detached state at 
any strength of adsorption e whereby for f = Jd no states with a particular value of m can be singled out as the most 
probable. Physically this means that for f ^ fo one expects very strong fluctuation of the number of contacs (which 
is our order parameter). In the following we show that this simplified consideration is indeed confirmed by the more 
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FIG. 1: (left) A snapshot from the MC simulation: iV = 128, h = 25.0, e = 4.0, ; (/) = 6.126. (right) Schematic representation 
of an adsorbed chain of length A'' which is pulled by the end segment off the surface with force /. While N = m monomers 
he on the plane, m monomers form the chain tail. The variation of F(m) with m G [l,Af] for super- f > fo and subcritical 
f < fo forces indicates two different minima (filled circles) of F at m = A'^ and m — for f > fo and f < fo, respectively. 

general adsorption model too. 



B. The Grand Canonical Ensemble approach to chain adsorption 

Starting with the conventional (i.e., force- free) case of polymer adsorption, we recall that an adsorbed chain is build 
up from loops, trains, and a free tail. One can treat statistically these basic structural units by means of the GCE 
approach [13, im where the lengths of the buildings blocks are not fixed but may rather fluctuate. The GCE-partition 
function is then given by an expansion over all possible lengths N, see Fig. ^j^, which can be considered and summed 
as a geometric series: 

In Eq. ([1]) 2 is the fugacity and U{z), V{z), and Q{z) denote the GCE partition functions of loops, trains and tails, 
respectively. The building block adjacent to the tethered chain end is allowed for by Vo{z) = 1-1- V{z). The partition 
function of the loops is defined as U{z) = {l^3z)"/n", where is the 3d connective constant and a is the 

exponent which governs surface loops statistics. It is well known that for an isolated loop a = 1 — 711 « 1.39 (l3| 
where 711 — —0.390. One can prove[14)| that a changes value, provided the excluded volume interactions between a loop 
and the rest of the chain are taken into account. The train GCE-partition function reads V{z) — {f^3wz)^/n^ 
with 1 — 7(i=2 ~ —0.343 whereby one assumes that each adsorbed segment gains an additional statistical weight 
w = exp(e). Eventually, the GCE partition function for the chain tail is defined by Q{z) = 1 + il^zz)"^ /n^ . For 

an isolated tail /3 = 1 — 71 « 0.32 where 71 — O.GSOfl^ but again the excluded volume interactions of a tail with the 
rest of the chain increase the value of (3. 

If one knows the GC partition function, Eq. ([T]), one can find the number of weighted configurations of a polymer 
chain, containing N segments (i.e., the canonical partition function of such chain), S^v, by taking the inverse Laplace 
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transform of 2(z). Using the generating function method [l5|, one finds that the main contribution to the coefficient 
Sjv at 2:^ is which is provided by the singularity at z* of S(z). There is a simple pole in Eq. (jT]) at 

z — z*, namely, when V{z*)U{z*) = 1. Thus one gets the free energy as = ksTN \n z* and the fraction of adsorbed 
monomers (which defines a convenient order parameter for the phase transition) is n = —dhiz*/dhiw. In terms of 
the so called polylog function, which is defined as z) = Yl^=i ^"/"■" JJJ and exists only for z <1, the equation 
for z* reads 

Af3Z*)$(A, = 1. (2) 

A nontrivial solution for z* in terms of w (or the adsorption energy e) appears at the critical adsorption point (CAP) 
w — Wc - see Fig. [2j) - where fji^z* = 1. For example, close to the CAP one may expand ^{a,n^z*) with respect to 
1 — ^^z* so that Wc is determined from C(Q;)<i>(l — jd=2, ^J'2Wc/ fJ-a) = 1 where ({x) denotes the Riemann ^-function. 
In the vicinity of the CAP the solution attains the form 



(a) 

Vo 



u 





u u 





(b) 


*(a,M,,z) 




'i '(X,ii,w'z ) 




* '(?t,H,w^z ) 




<1> '(^.lijW'z ) 









3 + 



FIG. 2: (a) Schematic representation of the series expansion, Eq. ([T}. (b) The intersection of the polylog functions ^{a, fi^z*) 
and 1/$(A, ^2'wz*) yields a solution of Eq. ([2} for the fugacity z* . For adsorption strength e < ec the corresponding Boltzmann 
weight ui^ = exp(e) is insufficient to provide an intersection point (the chain is desorbed) whereas for e > ec (for w^) a solution 
for z* exists (the chain is adsorbed). The CAP e = ec (i.e. for Wc) is marked by the first appearance of common point of 
intersection (full lines) at z* — l/ns- 



z*(w) « [1- A(w-u;,)i/("-i)]/i^' (3) 

where A is a constant. Then, for the average fraction of adsorbed monomers one obtains n oc (e — ec)^/'"~^-'~^. 
A comparison with the well known scaling relationship n oc (e — Ec)^^*^^^ where is the so called adsorption (or, 
crossover) exponent [l^ suggests that 

(j) = a-l (4) 



This result, derived first by Birshtein[ll|, is of principal importance. It shows that the exponent (f>, which describes 
polymer adsorption at criticality, is determined by the value of a which governs the polymer loop statistics! If loops 
are treated as isolated objects, then a = 1 — 711 w 1.39 so that = 0.39. In contrast, excluded volume interactions 
between a loop and the rest of the chain lead to an increase of a and 0, as shown below. 

From the expression for U{z), given above, and Eq. ([3]) we have Pioop ~ {p^z*y /l^^'^ « exp[— ci(e — ec)^/'^]/l^^'^ . 
This is valid only for e > Ec since a solution for Eq. ^ for subcritical values of the adhesive potential e does not exist. 
Nontheless, even in the subcritical region, e < Cc, the monomers occasionally touch the substrate, creating thus single 

I N-l 

loops at the expense of the tail length. The partition function of such a loop-tail configuration is Zi^t = ttt* jjfzjjw- 
On the other hand, the partition function of a tail conformation with no loops whatsoever (i.e., of a nonadsorbed 
tethered chain) is Zf = iV^^^^. Thus the probability P{^^p{l) to find a loop of length / next to a tail of length 

N — / can be estimated as P{^ap{l) = oc ^i+^^jy^^^^j for e < ec- In the vicinity of the CAP, e « Cc, the distribution 
will be given by an interpolation between the expressions above. Hence, the overall loop distribution becomes 

■jTT^exp [-ci(e-ec)^/'* , e > Ec 

^'loop(0 — \ FT? + ;i+*(jv-0f ' ^ — (5) 
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The same reasonings for a tail leads to the distribution 

e > Cc 

e = Cc (6) 
e < Cc 

In Eqs. ^ - ^ Ai, A2, Bi, B2 are constants. Close to the CAP these distributions are expected to attain a U - shaped 
form (with two maxima at Z « 1 and I « N), as predicted for a Gaussian chain by Gorbunov et al. [17]. 

For the average loop length L the GCE-partition function for loops yields L = zdU{z)/dz\z=z* = $(q; — 
l,^3Z*)/<I>(a,//3Z*). At the CAP, L diverges as L oc l/(e - eji/"^-!. 

The average tail length S is obtained as S" = zdQ{z) / dz\z=z' — ^{P — ^,^J■3Z*)/[l + $(/?,^3Z*)]. Again, using the 
polylog function, one can show that at tc the average tail length diverges as 5 oc l/(e — Cc)^^"^. 

C. The interaction of loops and the tail 

In the analytical expressions for the PDF of the different building units of a chain, Eqs. ([S])-® we didn't elaborate 
on the numerical values of the exponents a (that is, (j)) and (3, taking as an example those for non-interacting polymer 
chains. However, for a realistic self-avoiding chain one has to allow for the existence of excluded- volume interactions. 
To this end one may consider the number of configurations of a tethered chain in the vicinity of the CAP as an array 
of loops which end up with a tail. Using the approach of Kafri et al. Q along with Duplantier's [l2| graph theory 
of polymer networks, one may write the partition function Z for a chain with TV building blocks: TV — 1 loops and a 
tail[14j. Consider now a single loop of length M while the length of the rest of the chain is that is, M + K = N . 
In the limit of M > 1, A' > 1 (but with M/K < 1) one can show that Z - M''^^''^-i K^'m-i-^ where 
the surface exponent 7^ = 2-1- Af{v -\- 1) -\- ai + af and cti, trf are critical bulk and surface exponents [12]. The last 
result indicates that the effective loop exponent a becomes 

« = iH-i ~1H = ^ + ^ (7) 

Thus, 4> = a—\ = V = 0.588, in agreement with earlier Monte Carlo findings fT8|. One should emphasize, however, that 
the foregoing derivation is Mean-Field-like {Z appears as a product of loop- and rest-of-the-chain contributions) which 
overestimates the interactions and increases significantly the value of a, serving thus as an upper bound estimate. 
The value of a, therefore, is found to satisfy the inequality 1 — 711 < a <\ -\- i.e., depending on loop interactions, 
0.39 <cj)< 0.59. 

D. Taking the pulling force into account 

The GCE approach, described above, can now be employed to tackle the case of self-avoiding polymer chain 
adsorption in the presence of pulling force. Thus we extend the consideration of Gaussian chains by Gorbunov et al. 

As far as a force / is applied to the end-monomer of a tethered chain, one may choose two possible ways in which the 
chain detachment from the adsorbing surface can be carried out. One may fix / as an independent control parameter 
and study the variation of the height h of the end-monomer above the surface plane which corresponds to treatment 
within the constant force ensemble, herafter referred to as /-ensemble. Or, one might fix h and measure the force 
acting on the end- monomer at a given height, working thus in the constant height ensemble which we call in what 
follows the /i-ensemble. 



1. f -ensemble 

Under pulling force /, the tail GCE-partition function Q{z) in Eq. ^ has to be replaced by Q{z) = 1-1- 
y^^-AiiJ-^z)" ln^\ f d'^rPn{r) exp{fr±/T) where Pn(r) is the end-to-end distance distribution function for a self- 
avoiding chain [20] and fr± measures the work, spent to pull the chain end to height r± above the adsorbing surface. 
After some straightforward calculations Q{z) can be written as 



>exp[-ci(e-e,)V^ l] 

. lf(N-iy + 't' ■ 



Q(z) = 1 + ai /i-Ti/" $(1 - V, z^i^ exp{a2f'/n 



(8) 
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with the dimensionless force / = fa/ksT. The function Q{z) has a branch pomi singularity at z'^ = exp(— 02/^^"^), 
i.e., Q{z) ~ l/{z'^ — zY. One may, therefore, conclude that the total GCE-partition function ^(z) has two singularities 
on the real axis: the pole z\ related to the CAP, and the branch point z"^ , related to the pulling force. It is known 
(see, e.g.. Sec. 2.4.3. in [15|) that for 1 the main contributions to Sat come from the pole and the branch 

singular points, i.e., 

Ci{z*r'' + ^N''-'{z*r'' (9) 



Evidently, for large N only the smallest of these points matters. Note that z* depends on the adsorption energy e 
only (through w = exp(e)) whereas z"^ is controlled by the external force /. Therefore, in terms of the two control 
parameters, e and /, the equation z*{eD) = z'^^fo) defines the critical transition line between the adsorbed phase 



and the force-induced desorbed phase ■ 
it, / < /£), or above, / > /d, either z' 



Fig. [3l In the following this line will be called detachment line (DL). Below 
or z"^ , respectively, contribute to Sjv- The controU parameters, eo and /p, 
which satisfy this equation, denote detachment energy and detachment force, respectively. 

On the DL the system undergoes a first- order phase transition. The DL itself ends for /£> ^ in the CAP, Cc, 
where the transition becomes of second order, as is known for polymer adsorption without pulling. In the vicinity of 
the CAP the detachment force fo is predicted to vanish as fo'^i^ — ^cY^'^- 

This first order adsorption-desorption phase transition under pulling has a clear dichotomic nature (i.e., it follows 
an "either - or" scenario): in the thermodynamic limit N ^ 00 there is no phase coexistence! The configurations 
are divided into adsorbed and desorbed dichotomic classes. Metastable states are completely absent. Moreover, the 
mean loop length L remains finite upon DL crossing. In contrast, the average tail length S diverges close to the DL. 
Indeed, at f < fo the average tail length is given by = /i"^i/'^$(-z/, z*(w)/z#(/))/[l + ai$(l - i/, z*{w)/z* {/))]. 
At the DL, z* = " 
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FIG. 3: (a) Plot of the critical detachment force fo = fa/ksT against the surface potential e/ksT. Full and empty symbols 
denote MC and theoretical results. A double logarithmic plot of fo against e — with ec = 1.67 is shown in the inset, yielding 
a slope of 0.97 ± 0.02, in agreement with the prediction /o oc (e — ec)"^**. Shaded is shown the same phase diagram, derived by 
numeric solution of Eq. ^ along with z*{w) — z'^{f), which in dimensional / (right axis) against T (top axis) units appears 
reentrant, (b) The same phase diagram in units of detachment height hn and the distance from the CAP e — Ec. Dashed and 
solid lines denote theoretical predictions based on the Pincus, or Langevin force vs. elongation relationship while symbols show 
simulation data. 



2. h-ensemble 



In the constant height ensemble the way a chain tethered to a surface responds to stretching is described by the 
tail partition function Qn- The partition function of such chain with a fixed distance h of the chain end from the 
anchoring plane is given 

Qt.ii{N,h) = g^aPN{h) (10) 
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where again /3 = 1 — 71 and a is the bond length. The deformation of a polymer chain can be described within two 
models: the bead-spring (BS) model for flexible bonds and the freely jointed chain (FJ) model in which the bonds 
between monomers are considered rigid. In the BS model one can use for the expression [2l| 
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(11) 



where the exponent C w 0.8, and A is a normalization constant. The free energy of the tethered chain with a fixed 
distance h takes on the form Ftaii(-^, h) — ~T\nQtai\{N, h). By making use of Eqs. (fTO]) and (fTT|) . the expression for 
the force /jv, acting on the end- monomer when kept at distance h is given by 
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(12) 



One should note that at H/Rn ^ 1 we recover the well known Pincus deformation law: h oc aN{afN /kBTy-^'^^'^ ■ In 
this approximation the (dimensionless) elastic energy reads Ucx/ksT — —N{afpj/kBTy^'^. In result the corresponding 
tail free energy is given by 



Jtail 

ksT 



N 
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(13) 



Eq. indicates that there exists a height /iq = (C(l ~ v) / DY~'^ Rm over the surface where the force changes 
sign and becomes negative (that is, the surface repulsion dominates). According to Eq. (|12p the force diverges as 
/jv c)c —ksT/h upon further decrease of the distance h. 

It is well known [22] that the Pincus law, Eq. (fT2|) . describes the deformation behavior at intermediate force 
strength, l/A^" <C afN/ksT < 1. Direct Monte Carlo simulation results indicate that, depending on the model, 
deviations from Pincus law emerge at H/Rn > 3 (bead-spring off-lattice model) [l^], or H/Rn > 6 (Bond Fluctuation 
Model) . In such "overstretched" regime (when the chain is stretched close to its contour length) one should take 
into account that the chain bonds cannot expand indefinitely. This case could be treated within the simple FJ model 
[2^ 1 where the bond length a is fixed. In this model the force - elongation relationship is given by 



knT 
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(14) 



where C ^ denotes the inverse Langevin function C{x) = coth(2;) — 1/x. The corresponding free energy of the tail for 
the FJ model reads 



^tail 

k^ 



= -NQ 



N 



T 



(15) 



where we have used the notation G{x) = xC{x) = a;coth(x) — 1. One should emphasize that the force Jn stays 
constant in the course of the pulling process (i.e., as long as one monomer, at least, is adsorbed on the surface), 
thus /at corresponds to a plateau on the elongation curve f — h. An adsorbed monomer has a chemical potential, 
Mads = In 2*, which should be equal in equilibrium to the chemical potential of a desorbed monomer in the tail, 
Mdes = d{Ftaii/T)/dN. Thus the condition /Xads = A^dos leads to the following "plateau law" relationship 
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(16) 



with Q 1 being the inverse of the g function. Close to the critical point ec the plateau force /p — > 0. Indeed, taking 
into account that in the vicinity of the critical point ln[/i3Z*(e)] oc — (e — ec)^/'^ and G~^{x) ~ (Sa;)^/^ we conclude 
that /p oc (e - Cc)''/'^ for the BS model and /p oc (e - e^Y^^'t' for the FJ model. If the number of tail monomers is 
denoted by M, then the one can write 14] n ~ —{l/TN)dFg^^s/de, where i^ads is the free energy of the adsorbed 
portion of the chain given as Fads = kBT[N — M(/i, e)] lnz*(e). From Eq. ^TE\i one can easily obtain M for h ^ Rg 
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so that in result one gets 
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(17) 



where Ci, C2, C3 are constants of the order of unity. The derivatives C'{x) = — l/[sinh(x)]^ and G'{x) = 

£{x) + xC'{x) = coth(x) — x/[sinh(x)]^ and ci = (1 — i^)ci. 

As one can see from Eq. ()17p . the order parameter decreases linearly and steadily (no jump!) with growing h/N . 



E. Reentrant phase behavior 

Recently, it has been realized (25l | that the DL, force Jd versus temperature T, when represented in units with 
dimension, goes (at a relatively low temperature) through a maximum, i.e., the desorption transition shows reentrant 
behavior! Such behavior has been predicted earlier [2^ [23, [IM] in a different context, namely, of DNA-unzipping, and 
also in the coil- hairpin transition[29|. 

One can readily see that this result follows directly from our theory. Indeed, the solution of Eq. ^ at large values 
of e (that is, at low temperature) can be written as z* k, e~'^ / 113 so that the DL, z* — z^^, in terms of dimensionless 
parameters is monotonous, fjj oc [eo — Iii(/i3//i2)]''- Note, however, that the same DL, if represented in terms of 
the dimensional control parameters, force fo versus temperature Tjj (with a fixed energy Eq), shows a nonmonotonic 
behavior Jd = kBToleo/TD — \n{fj,3 / ii2)Y / a - Fig. [3l as found earher for DNA-unzipping This curve has a 

maximum at a temperature given by ksTj^"'^ = (1 — i/)eo/ ln(/X3//i2). At very low T, however, the expression for 
P„(r) [ipj predicts divergent chain deformation [2a], i.e., it becomes unphysical. One can readily show that in this 
case the correct behavior is given by fa — Eq + kBT\n.{iJ,^/ ijl2)- 



III. MONTE CARLO SIMULATION MODEL 

We use a coarse grained off-lattice bead-spring model [s^ which has proved rather efhcient in a number of polymers 
studies so far. The system consists of a single polymer chain tethered at one end to a flat impenetrable structureless 
surface - Fig. [TJ The surface interaction is described by a square well potential. 



e, z < rc 
0, z > r. 



(18) 



The strength e is varied from 1.0 to 7.0 while the interaction range rc — 0.125. The effective bonded interaction is 
described by the FENE (finitely extensible nonlinear elastic) potential: 



Ufene = -K{1 - a^ln 



I — a 



(19) 



with K ^ 20, Iraax 

Morse potential: 



1, a = 0.7, l„ 



0.4. The nonbonded interactions between monomers are described by the 



UM{r) 

CM 



exp(-2a(r - r,„i„)) - 2exp(-a(r - r-mm)) 



(20) 



with a ~ 24, rmin — 0.8, eM/ksT — 1. In few cases, needed to clarify the nature of the polymer chain resistance 
to stretching, we have taken the nonbonded interactions between monomers as purely repulsive by shifting the Morse 
potential upward by ejvf and removing its attractive branch, Vj^ir) — for r > r^in. 

We employ periodic boundary conditions in the x — y directions and impenetrable walls in the z direction. The 
lengths of the studied polymer chains are typically 64, and 128. The size of the simulation box was chosen appropriately 
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to the chain length, so for example, for a chain length of 128, the box size was 256 x 256 x 256 . All simulation results 
have been averaged over about 2000 measurements. 

The standard Metropolis algorithm was employed to govern the moves with self avoidance automatically incorpo- 
rated in the potentials. In each Monte Carlo update, a monomer was chosen at random and a random displacement 
attempted with Ax, Ay, Az chosen uniformly from the interval —0.5 < Ax, Ay, Az < 0.5. If the last monomer 
was displaced in z direction, there was an energy cost of — /Az due to the pulling force. The transition probability 
for the attempted move was calculated from the change AU of the potential energies before and after the move was 
performed as = exp{— AU / kBT) . As in a standard Metropolis algorithm, the attempted move was accepted, if W 
exceeds a random number uniformly distributed in the interval [0, 1]. 

As a rule, the polymer chains have been originally equilibrated in the MC method for a period of about 5 x 10^ 
MCS after which typically 500 measurement runs were performed, each of length 2 x 10^ MCS. The equilibration 
period and the length of the run were chosen according to the chain length and the values provided here are for the 
longest chain length. 



IV. COMPARISON OF SIMULATION DATA WITH THEORETICAL PREDICTIONS 



We have investigated the force induced desorption of a polymer performing MC simulations in the /-ensemble 
and in the /i-ensemble. As an order parameter for the desorption transition we use the fraction of monomers n in 
contact with the sticky surface. Below we present few typical quantities of interest which manifest the good agreement 
between theoretical predictions and simulation results. Another important point is the observed qualitative difference 
between the /— and ft,— ensembles in the behavior of some basic properties like the order parameter of the phase 
transition. Fig. 3^ shows the variation of the order parameter n with changing adhesive potential e in the /-ensemble 
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FIG. 4: (a) The 'order parameter', n, against the surface potential, e, for various pulling forces /. The chain has length A'^=128. 
(b) Plot of n vs. / for several surface potentials e. 



at fixed pulling force whereas FigHIa depicts n vs. force fa/T for various e. The abrupt change of the order parameter 
is in close agreement with our theoretical prediction. Indeed, from Fig. [4] one can readily verify that the polymer 
detachment transition is of first order. 

However, the order parameter variation in the equivalent /i-ensemble looks very different. In Fig. [5^, [SJa, we show 



the change in n with h and in the insets the variation of the fraction of adsorbed segments with adsorption strengths 
e for several fixed heights 20 < /i < 50 of the N — 128 chain. It is evident that, apart from the rounding of the 
MC data for n at n ^ 0, which is less pronounced for N — 128 than for N — 64, one finds very good agreement 
between the behavior, predicted by Eq. p?]) . and the simulation results. Comparing Figs, d] and [5] one realizes the 
striking difference between the order parameter behavior in the /— and ft— ensembles. However, if the height ft on 
the X-axis of Fig. [5k is expressed in terms of the corresponding average force (/), one recovers again a jump in the 
order parameter n [8|. 

The peculiar nature of the desorption transition under pulling becomes more evident when one plots the PDF of 
the order parameter in both statistical ensembles. In the presence of a pulling force one observes a remarkable feature 
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FIG. 5: (a) Order parameter n variation with changing height h/aN of the fixed chain-end for polymers of length A*' = 64, 128 
and different adsorption strength e/ksT. (b) Variation of n with e/ksT for different fixed positions of the chain-end h/aN. 



of the order parameter probability distribution - Fig. [5^: - an absence of two peaks in the vicinity of the transition 



0.07 




FIG. 6: Order parameter distribution of a polymer with A'^ = 128 and ads adsorption potential e = 3.0: (a) in the /—ensemble 
for several values of /. The critical detachment force is fn ~ 1.85 ± 0.01. (b) in the /i— ensemble for several fixed heights h. 
The inset displays the distribution W{n) at the critical height of detachment ho ~ 54.3. 

force although bimodality is customary in first-order phase transition. Immediately at Ju the distribution W{n) 
is flat, indicating huge fluctuations of n so that any value of the number of contacts is equally probable. This lack of 
bimodality in the W{n) manifests the dichotomic nature of the desorption transition which rules out phase coexistence. 
In contrast, in the ft,— ensemble, Fig.[6}3, one observes an entirely different shape of W{n) with only slight deviations 
(an appearance of non-zero third moment of the distribution) from Gaussianity in the vicinity of ho ■ The fluctuations 
of n, according to the half-width of W{n), remain flnite and almost unchanged for all values of h. 

Eventually, we show in Fig. [7^ the typical plateau observed in the average pulling force when the polymer detachment 
is effected in the /i-ensemble. Within a large interval of height variation the mean force, exerted on the end monomer, 
remains constant as observed in laboratory experiments. A rapid growth in the magnitude of this force sets in after 
the plateau, as soon as the bonds rather that the confortmation of the polymer are stretched upon further elongation. 
The stronger the adsorption, e, the larger the force fo required to remove the chain from the substrate. 

In addition to the force due to bonded interactions, however, one can see a small contribution from the non-bonded 
(attractive) interactions between the chain segments. This contribution is not allowed for by the GCE theory and, 
therefore, a test with the theoretical preedictions should exclude it. If the attractive branch of the Morse potential is 
removed, leaving the self-excluded repulsive branch only, this contribution almost vanishes - Fig. [7^ (inset). 

The elongation vs. force relationship, predicted by Eq. ((H]), is tested in Fig.[7lD for chains in which only non-bonded 
repulsion between segments exists. For small and intermediate extensions h where / is not too large the agreement 
with Pincus law is found to be perfect although it deteriorates for larger /, as expected. In the latter region one may 
show that a very good agreement between theory and computer experiment is provided by the FJ model - Eq. (jl4p . 
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FIG. 7: (a) Variation of the two components to the total force, exerted by the chain on the end-monomer which is fixed 
at (dimensionless) height h/aN for diff'erent adsorption potentials 2.0 < e/ksT < 5.0 and bonding (FENE) interactions (full 
symbols) as well as non-bonding (Morse) interactions (empty symbols). In the inset the same is shown for a neutral plane 
e = 0.0 and purely repulsive monomers (triangles) and for the usual Morse potential (circles), (b) Variation of the total applied 
force / with growing height of the end monomer in terms of Pincus reduced variables, faN'^ /ksT versus h/aN" , for a polymer 
with purely repulsive nonbonded forces for A'' = 64, 128. 



From Fig. [TJs one can also see that / goes through zero before the height has become zero, that is, no force is felt 
when the chain end is kept at this particular height. Further decrease of h leads to change of sign of /, indicating the 
entropic repulsion of the polymer coil from the solid surface. 



V. SUMMARY 



In conclusion, we have shown that a full description of the force-induced desorption of a self-avoiding polymer chain 
can be achieved by means of the GCE approach, yielding the average size and probability distribution functions of 
all basic structural units of partially adsorbed polymer as well as their variation with changing force or strength of 
adhesion. All these predictions appear in good agreement with our MC simulation results. 

The polymer detachment transition under pulling is found to be of first order whereby due to its dichotomic nature 
phase coexistence is impossible. This absence of binodal states makes the polymer desorption under pulling a rather 
unusual in comparison to conventional first-order phase transformation. 

The critical line of desorption, while growing steadily when plotted in dimensionless units of detachment force against 
surface potential, appears "reentrant" in absolute units of force against temperature. Thus, at very low temperature 
the polymer is expected to be desorbed, with the growing T it may adsorb, and at even higher temperature - desorb 
again from the surface. 

One finds that the crossover exponent, (/), governing polymer adsorption at critical! ty, whose exact value has been 
controversial for a long time, depends essentially on interactions between different loops so that <j) may only vary 
within the hmits 0.39 < (j) < 0.59. 

A point of more general importance for the statistical mechanics in general and theory of phase transitions in 
particular is the issue of ensemble equivalence. The latter implies an identity of the equation of state, regardless of 
which statistical ensemble has been employed, whereas the fluctuations within the different ensembles may be entirely 
different Q. For finite polymer lengths, however, differences in the equation of state may also be visible. As far as 
in practice one deals with finite polymer chains in laboratory experiments, this difference is expected to be clearly 
manifested in cases of practical concern. 
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